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ABSTRACT 

Spiral structure is the most distinctive feature of disk galaxies and yet debate persists about which theory of 
spiral structure is correct. Many versions of the density wave theory demand that the pitch angle be uniquely 
determined by the distribution of mass in the bulge and disk of the galaxy. We present evidence that the tangent 
of the pitch angle of logarithmic spiral arms in disk galaxies correlates strongly with the density of neutral 
atomic hydrogen in the disk and with the central stellar bulge mass of the galaxy. These three quantities, when 
plotted against each other, form a planar relationship that we argue should be fundamental to our understanding 
of spiral structure in disk galaxies. We further argue that any successful theory of spiral structure must be able 
to explain this relationship. 

Subject headings: galaxies: evolution — galaxies: spiral — galaxies: structure — galaxies: fundamental 
parameters — galaxies: kinematics and dynamics 


1. INTRODUCTION 


Spiral structure is a commonplace and visually strik¬ 
ing feature of many galaxies and yet there is still dis¬ 
agreement as to the correct theory that explains its ori¬ 
gin after decades of debate. The first well known the¬ 
ory dLin & Shul H%4h is that density waves propagating 
through the disk of the galaxy are the responsible agent. 
The density wave theory for spiral modes, described in de¬ 
tail by iBertin & Linl (119961) . calls for a long-lived, quasi¬ 
steady global spiral pattern. Others call for more tran- 
sient spiral patterns, whether from swing-amplified n oise 
dGoldreich & Lvnden-Belliri965l : I Julian & Toomre|| 19661) . re¬ 


current cycles of groove mode s (ISellwood 2000), or su¬ 
perposed transient instabilities (ISellwood & Carlberdl2014 ). 
Other theories have also been proposed, with one in particu - 
lar, the manifold theory (lAthanassoula et al.ll2009bllaL l2010t) . 


rejecting the density wave concept altogether in favor of an 
explanation involving stars in chaotic highly eccentric orbits. 

The density wave theory, as originally articulated by 
lLin&Shul(ff966l) . had a very specific prediction for the pitch 
angle of the spiral pattern produced by the waves. They cal¬ 
culated the pitch angle to be a ratio of the density of material 
in the galaxy’s disk to a certain quantity made up of the fre- 
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quencies of orbital motions in the disks, 

2irmG(<r 0 + F cr*) 
tan|F| = «(V- (a ,-,„iW 


( 1 ) 


where P is the logarithmic spiral arm pitch angle, m is the 
number of spiral arms in the pattern, G is the gravitational 
constant, O is the angular frequency of orbits of particles in 
the disk, R is the galactocentric radius, k is the epicyclic fre¬ 
quency of the same particles, and ui is the frequency asso¬ 
ciated with the perturbation that excites the density waves. 
Note that as long as this perturbation is some form of self¬ 
excitation arising from within the disk itself, it follows that all 
of the terms in the denominator should depend on the mass 
of the central bulge of the galaxy. The simplest case of a 
dominant bulge (approximated as a point mass) would im¬ 
ply, for instance, f l ex \J~M and n oc \/M, with M being the 
central mass. The numerator depends on the density of gas 
in the disk, a 0 , and the stellar disk density, cr„, with a factor 
F, called the reduction factor, which underweights the stel¬ 
lar density (compared to the gas density), since it is primarily 
within the gas that the density wave propagates. In this Letter 
we will present some evidence that F << 1. 

Focusing on the masses and densities involved in this rela¬ 
tion, we find that 


tan|R| oc 


<J 0 + Fa* 

M n 


( 2 ) 


where M a is the mass of the galactic bulge, or else the total 
mass interior to the radius in question. 

This formula is known to work very well in the case of 
spiral density waves in Saturn’s rings dShul 119841) . Bulge- 
dominated galaxies are not too distant from the Saturnian sit¬ 
uation of a small dense core with negligible mass in the disk, 
though disk-dominated galaxies are obviously far more com¬ 
plex. Generally speaking, the density wave theory predicts 
that the pitch angle of the spiral arms in galaxies does de¬ 
pend on the radial distributi on of matter in the galax y, and 
experimental studies concur (iSeigar et al.|[200 61 120141) . Re¬ 
sults of this type broadly agree with lLin & Shul (11966 1) that a 
thin (dense) disk and massive (small) central bulge should re¬ 
sult in a tight (loose) spiral. This is not surprising, since we 
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would expect a standing wave pattern (such as in a vibrating 
string) to depend on the ratio of a restoring force or tension 
(in this case the central mass, or at least the mass inside a 
given radius R ) to the density of the medium (in this case the 
density of gas in the disk at radius R). Although the precise 
nature of the relation between these three quantities can be 
expected to vary between galaxies of different types (bulge- 
dominated versus disk-dominated, for instance), nevertheless 
we show in this Letter that the three quantities, spiral arm 
pitch angle, central bulge mass, and gas density in the disk, do 
strongly correlate to form a fundamental plane that may play 
a similar role in tying together gross features of disk galaxies 
to that played by the fundamental plane of ellip tical galaxies 
(iDiorgovski & Davisll 1987k iDressler et a 1.11 1 9871) . 

We take as our sampl e disk galaxies from the DiskMass 
Survey (DMS: iBershadv et al.112010 ). which is ideal for our 
purposes since it deals with the disk densities of a sample 
of face-on galaxies and includes measurement^of the central 
bulge mass. Using the technique of iDavis et al.l (120121) . we 
measure the pitch angle for these galaxies and find that our 
sample of 24 galaxies, when plotted in a volume defined by 
these three quantities, delineates a plane with very low scat¬ 
ter. There is only a 0.0047% chance that this plane could have 
been formed by statistical accident. 

The plane satisfies a number of requirements, which one 
would expect of a useful fundamental plane. The plane is 
steeply inclined across the volume formed by the three re¬ 
lated quantities. In other words, it is not merely a relation 
between two of the three quantities, with the third essentially 
irrelevant. The galaxies are distributed quite widely and fairly 
uniformly across the plane. There is no particular evidence of 
a favored curvilinear relation on the face of the plane. Finally, 
and most importantly, the plane is oriented as one would ex¬ 
pect on the basis of the density wave theory. A large bulge 
mass and a rarified disk produces the tightest spirals. A small 
bulge and a dense disk produce the loosest spirals. We submit 
that any successful theory of galactic spiral structure must be 
able to explain this result. 

One final point is worthy of note. The DMS measured 
not only the density of atomic hydrogen in the disk of each 
galaxy (the quantity used in our relation), but also the density 
of molecular hydrogen and the dynamical disk mass density 
(the total density in the disk). Our results suggest that it is 
the gas density, not the total density in the disk that matters 
for spiral density waves. This suggests an apparent decou¬ 
pling between the stars and gas. The fact that the density of 
atomic hydrogen fits noticeably better than that for molecular 
hydrogen may simply be due to the fact that it is a much more 
reliable measurement, since molecular hydrogen is estimated 
indire ctly from observations of other molecules, not hydrogen 
itself (I Westfall et al .1120lit iMartinsson et al 12013 all . 

2. DATA AND ANALYSIS 

The DMS PPak Sample (IMartinsson et al.l[2013b ) consists 
of 30 nearly face-on galaxies whose disk densities have pre¬ 
viously been closely studied. However, four of these galaxies 
do not have central stellar bulge masses available and three 
provide them only as upper limits, so they are excluded from 
our sample. In addition to the 23 remaining DMS galaxies, we 
also include our own Galaxy, the Milky Way, in our sample of 
galaxies. 

In addition to the method described in IDavis et al.l (120121) . 
which utilizes a two-dimensional fast Fourier transform soft¬ 
ware called 2DFFT, we also measured pitch angles for all of 


the sample galaxies using new software called Spirality. Spi- 
rality (Shields et al. 2015, in preparation) is a novel method 
for measuring spiral arm pitch angle by fitting galaxy images 
to spiral coordinate systems (templates) of known pitch. For 
a given pitch angle template, the mean pixel value is found 
along each of typically 1000 spiral axes. The fitting function, 
which shows a local maximum at the best-fit pitch angle, is the 
variance of these means. In other words, we choose the pitch 
angle that exhibits the greatest contrast between the mean lu¬ 
minosity along the spiral axes. The presumption is that where 
the pitch angle of the spiral axes is equal to the pitch angle of 
the galaxy’s spiral arms there will be some axes that fall pre¬ 
cisely along the true spiral arms (and thus are much brighter 
in the mean) and some that never coincide with the true spiral 
arms (and thus are, on average, dim). Where the pitch of the 
axes is not equal to the pitch of the spiral arms, each axis will 
cross the true spiral arms a roughly equal number of times, 
making the mean brightness along each axis roughly equal. 
Error bars are found by varying the inner radius of the mea¬ 
surement annulus and finding the standard deviation of the 
best-fit pitch angles. 

The two techniques yield measurements that agree within 
the error bars in almost all cases. As a final and important 
test, we visually inspected each galaxy, comparing them to 
overlays of synthetic spirals on transparency paper, in order to 
confirm the measured pitch angle. Our overlay transparencies 
showed spirals of different sizes and different pitch angles in 
steps of 5° from 5° to 85°. We were therefore able to visually 
confirm the pitch angle to within 5°. We were satisfied in all 
cases that the measured pitch angle of 2DFFT was reliable 
and strongly supported by the combination of Spirality and 
visual inspection. For the sake of consistency, we chose to 
use only the results of 2DFFT in this Letter. The pitch angles, 
P, given in Table Q] are the results of the 2DFFT routine. 

The images used were obtained from the NASA/IPAC Ex- 
tragalactic Database^] and/or from the pODI (partial One 
Degree Imager) camera on the WIYN 3.5 m telescope. The 
WIYN images were all acquired as 120 s exposures, cali¬ 
brated using OuickReduce 1.0 from the ODI Pipeline, Portal, 
and Archive FI and processed using a five-point dither pattern 
for each g alaxy and subse quently stacking the images using 
SWarp dBertin et al.112002) . Additionally, KPNO 2.1 m imag¬ 
ing for UGC 463, 1529, 1908, 4036, and 11318 were mea¬ 
sured to confirm previous pitch angle measurements. Unless 
otherwise specified (Milky Way data have been determined 
in very different ways than other galaxies), all data for stellar 
galactic bulge masses (M*“ ,g< ’) and maximum neutral atomic 
hydrogen (Hi) gas ma ss surface densities (£™(“ ) come from 
iMartinssonl (1201 ll) and lMartinsson et alJ (120 1 3al) : see TableQ] 
For the determination of the stellar bulge masses, the AT-band 
light profile was decomposed into a central Sersic component 
(conv olved with a seeing disk) and a number of exponential 
disks (IMartinsson et al.ll2013bb . The bulge masses were de¬ 
termined using the integral of the light from the central Ser¬ 
sic component and the mass-to-light ratio ( M/L ) derived from 
the dis k using vertical velocity dispersions (IMartinsson et al.l 
1201 3al) . Ga s densities were d etermined from 21 cm line mea¬ 
surements (lMartinssonll201 ll) . 


10 http ://ned. ipac .caltech.edu/ 

11 http:/ /portal. odi .iu.edu 
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Table 1 

Sample 


Galaxy Name 
(1) 

Type 

(2) 

Band 

(3) 

Image Source m 
(4) (5) 

tan |P| 

(6) 

\og(M b : ,,se /M q) 
(V) 

£™7(M 0 pc“ 2 ) 

(8) 

Excluded 

(9) 

Milky Way 

SBc 

21 cm 

1 

4 

0.414 

± 0.051 a 

9.95 ± 0.03 b 

4.98 

±0.53 c,d 


UGC 448 

SABc 

r 

2 

4 

0.327 

±0.033 

9.76t 23 

4.58 

±0.46 


UGC 463 

SABc 

B 

3 

3 

0.412 

±0.066 

9.35(ffj 

6.18 

±0.66 


UGC 1081 

SBc 

r 

2 

2 

0.452 

±0.064 

°' 81 -0.24 

6.25 

±0.62 


UGC 1087 

Sc 

r 

2 

2 

0.188 

±0.039 

8-64^i| 

4.42 

±0.46 


UGC 1529 

Sc 

645.0 nm e 

4 

3 

0.490 

±0.096 

8.987 2 ° 

6.59 ±0.66 


UGC 1635 

Sbc 

r 

2 

3 

0.209 

±0.014 

8 74^ 

5 ■'^-0.29 

2.60 

±0.32 


UGC 1862 

SABcd 1 

r 

2 

2 

0.444 

± 0.074 


9.14 

±0.91 

/ 

UGC 1908 

SBc 2 

645.0 nm e 

4 

3 

0.376 

±0.069 

9.68t 23 

4.62 

±0.46 


UGC 3091 

SABd 

i 

2 

2 

0.555 

±0.092 


5.59 

±0.56 

/ 

UGC 3140 

Sc 

r 

2 

3 

0.290 

±0.090 

9.65t>l 

4.87 

±0.54 


UGC 3701 

Scd 

r 

2 

2 

0.276 

±0.090 

8.69(9!? 

5.55 

±0.57 


UGC 3997 
UGC 4036 

Im 

SAB be 

g 

645.0 nm e 

5 

4 

2 

2 

0.185 

0.268 

±0.048 

±0.021 

8 53+^ 

8.92^ 

5.01 

5.20 

±0.54 

±0.56 


UGC 4107 
UGC 4256 

Sc 

SABc 

g 

g 

5 

5 

2 

2 

0.371 

0.555 

±0.041 

±0.099 

8 65+&S 

Q 7Q+0.36 
^• z ^-9.29 

5.42 

9.75 

±0.54 

±0.98 

/ 

UGC 4368 

Scd 

g 

5 

2 

0.439 

±0.043 

9-2KP 

5.95 

±0.66 


UGC 4380 

UGC 4458 
UGC 4555 

Scd 

Sa 

SABbc 

g 

g 

g 

5 

5 

5 

3 

1 

2 

0.430 

0.243 

0.213 

±0.095 

±0.056 

±0.017 

8-86 

10 67 +(a -° 

au - d/ -Q.39 

8 96 +0 ^° 

°- yo -0.39 

4.08 

3.28 

4.58 

±0.41 

±0.53 

±0.47 


UGC 4622 

Scd 

g 

5 

4 

0.401 

±0.099 

9.89tii 

3.50 

±0.38 


UGC 6903 

SBcd 

g 

5 

2 

0.283 

±0.041 

8' 03 -q il 

4.94 

±0.59 


UGC 6918 

SABb 3 

F606W 

6 

3 

0.306 

± 0.044 

8.04;}$ 

7.04 

±0.72 

/ 

UGC 7244 

SBcd 

g 

5 

2 

0.627 

±0.105 


5.53 

±0.60 

/ 

UGC 7917 
UGC 8196 
UGC 9177 

SBbc 

Sb 

Scd 

g 

g 

g 

5 

5 

5 

3 

5 

2 

0.278 ±0.025 
0.144 ±0.009 
0.256 ±0.035 

10 oi +00 - 34 

iu.ui_, 0 ni 

10 73+ 016 

AU -' -0.26 

Q cc+0.2-4 

2.70 

2.74 

3.92 

±0.28 

±0.28 

±0.42 

/ 

UGC 9837 

SABc 

g 

5 

6 

0.482 

±0.061 

O T2C+0.17 
°.jj_029 

7.95 

±0.80 


UGC 9965 

Sc 

g 

5 

3 

0.237 

±0.037 


5.63 

±0.58 

/ 

UGC 11318 

SBbc 

645.0 nm e 

4 

3 

0.569 

±0.101 

9-69t| 

6.51 

±0.67 


UGC 12391 

SABc 

r 

2 

4 

0.235 

±0.091 

8-98}” 

4.90 

±0.49 



Note. — Columns: (1) Galaxy name. (2) Hubble morphological type from either the UGC dNilsonl 1973D or RC3 dde Vaucouleurs et alj 
fTWIh catalogs. Notes on morphologies: 1 = peculiar, 2 = starburst, and 3 = AGN. (3) Filter waveband/wavelength used for pitch angle 
calculation. (4) Telescope/literature source of imaging used for pitch angle calculation. (5) Harmonic mode (number of spiral arms). (6) 
Tangent of the pitch angle of the galactic logarithmic spiral arms. (7) Base 10 logarithm of the stellar bulge mass of the galaxy, in solar 
masses. (8) Maximum surface density in the galactic Hi gas, in solar masses per square pc. (9) Indicatio n of galaxies that are excluded in 
fittings due to missing measurements or measurements that are merely upper limits. Image Sources: (l) ILevine et all (120061) : (2) WIYN 
3.5 m pODI; (3) JKT 1.0 m; (4) Palomar 48 inch Schmidt; (5) SDSS; (6) HST. 
a Levine et aD <20061) . 
b lMcMillanll201lh . 

c No error estimates were provided by its reference so we have assigned the mean error of the included sample, ±0.53 Mqpc~~. 
d Calculated using Equation (2) from lFerrierel <20011 1. 
e IllaJ emulsion. 


3. RESULTS 

We find a best linear fit for Equation (0 from the included 
data sample of 24 galaxies of 


tan |P| = (0.375 ±0.092) 


S^V(M 0/ 7c- 2 ) 

log(M^7M 0 ) 


±(0.127 ±0.049). 


(3) 

The root mean squared error ( E nns ) is equal to 0.0909 (a resid¬ 
ual scatter of 31.2% per galaxy on average), with R 2 = 0.344, 
and a /t-value equal to 2.59 x 10~ 3 for Equation (0. A plot 
of this linear fit, along with the included data sample, is given 
in Figure [I] 

The formula describing the fundamental plane for spiral 


galaxies from the sample is as follows: 

Y^nax 

-——y = (5.70 ± 1.40) tan IP I- 

Mqpc 2 

(0.677±0.199)log(M7 s 7M o )±(9.29± 1.96). (4) 

E nn s = 0.770 Mqpc~ 2 (a residual scatter of 16.7% per galaxy 
on averagJT with R 2 = 0.613, and a p-value = 4.71 x 10 -5 
for Equation 0. It is interesting to note that the addition of 
the extra dimension cuts the residual scatter approximately in 
half. A three-dimensional plot of this plan e, along with the 
included data sample, is given in Figure 00. 

12 Compare this to the fundamental plane for elli ptical galaxies, which has 
a residual scatter of ~ 20% per galaxy on average iK ormendv & D iorgovskil 
09891) . 

13 A 3D animated gif of this figure can be accessed at 
http://dafix.uark.edu/~ben/movie.gif 
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Figure 1. Two-dimensional plot of the linear fit defined by the multivariate 
normally distributed sampling fit of Equation (3) depicted by the solid black 
line, along with the plotted points of the 24 galaxy member data set. The 
Milky Way is depicted distinctly in green. The axes [x,y] depict [tan|P|, 
(jy^ax / (Mq pc ~ 2 )) / \og(M^ ulge /M q )], respectively. 

The errors presented in Equations (0 and 0 have been de¬ 
termined by sampling individual data points from multivariate 
normal distributions using the following algorithm. 

1. For each measurement, draw a new measurement based 
on multivariate normal distributions with the mean and 
variance of each variable for all 24 galaxies. 

2. Fit linear (or planar) best-fit coefficients to the 24 
points. 

3. Repeat steps 1 and 2 10 6 times, saving the fitted coeffi¬ 
cients after each step. 

4. Use the distribution of the coefficients based on these 
10 6 fits to define the best-fitting (median) coefficients 
and their error (the la confidence interval of the distri¬ 
bution). 

The orientation of the fundamental plane illustrated in Fig- 
ure[2]is exactly as one would expect on the basis of the spiral 
density wave theory. According to Equation 0, the pitch an¬ 
gle is minimized (tightest winding) when the Hi mass surface 
density is low and the central mass is high. Alternatively, the 
pitch angle is maximized (loosest winding) when the H t mass 
surface density is high and the central mass is low. This be¬ 
havior is illustrated in both figures. Note in the middle panel 
of Figure [2] how the plane slopes from the lower left front 
(low pitch angle, low Hi mass surface density, and high cen¬ 
tral mass) corner of the cube to the upper right back (high 
pitch angle, high H : mass surface density, and low central 
mass) corner of the cube. Furthermore, this indicates that the 
shape of the plane is strongly correlated to all three variables 
(the individual variable p-values of the intercept, tan |/ J |, and 

log(M*'^7M 0 ) are 9.97 x 1(T 5 , 8.03 x 10“ 5 , and 9.41 x 
10 -4 , respectively for Equation 0). 

4. DISCUSSION 

If one favors the standing wave picture of spiral structure, 
our result is not unexpected. In analogy with standing waves 
on a string, we would expect the wavelength of the resulting 
pattern to be strongly determined by the tension or restoring 
force (in this case the central gravitating mass) and the den¬ 
sity of the medium (in this case the gas in the disk). It is 


worth noting that our case is probably analogous to a string 
with non-uniform density, since the gas density generally falls 
off with increasing radius in a galactic dislo Additionally, 
it seems reasonable that the gravitational restoring force in¬ 
creases with increasing radius, since there will be more mass 
inside the given radius. Both effects would tend to cause the 
pitch angle to tighten with increasing radius,jtnd this effect is 
often seen in spiral galaxies (see, e.g., iDavis et a02012l) . 

It has been often proposed in the past that different mech¬ 
anisms may explain spiral structure in different galaxies; for 
instance, the mechanism that produces grand design spirals 
ma y differ from the one that produces flocculent spirals (e.g., 
iD’Onghia et al.l20l3l) . Furthermore, it has been demonstrated 
that galaxies that appear to have grand design structure in in¬ 
frared light can appear flocculent in blue images , sug gesting 
that stellar and gaseous disks are decoupled dGrosbol & Patsisl 
119981) . It is striking that the sample used in this study contains 
quite a few flocculent or multi-armed patterns, which are not 
grand design. The existence of a very low scatter planar cor¬ 
relation for all of these galaxies is thus very significant and 
implies that different galactic morphologies all adhere to the 
same imposed mechanism of density wave theory. It is true 
that four of the galaxies have a noticeably greater scatter than 
the others, and further study with larger samples might yet 
support the existence of two kinds of spiral structure. It is 
worth noting that two of these four galaxies represent the ex¬ 
tremes of gas density for the sample, one having clearly the 
highest gas density in its disk, another clearly the lowest. 

In recent years, there has been some discussion that spi¬ 
ral arms may be quite transient, persisting for only one o r 
two revolutions of the disk galaxy (iToomre & Kalnai sill 9911) . 
In recent years, there have been attempts to show theo- 
retically that more long-lasting spiral patterns ar e possible 
(ISellwood & Carlbergll984[D’Onghia et al.l2013i) . This Let¬ 
ter suggests that even if spiral patterns are transient, some res¬ 
onant mechanism compels the pattern, when it reforms, to re¬ 
sume something close to its previous pitch angle. 

The relation discovered here might be useful as a tool in 
the study of disk galaxies. One of the three quantities, disk 
density, is relatively difficult to measure. The relation found 
here could be used to measure it indirectly from the other two 
quantities (pitch angle and central bulge mass), which would 
be easier to measure. In addition, the existence of the three- 
way correlation may enable more careful studies of the im¬ 
portant relation between pitch angle and central mass, which 
is itself a very us eful marker for quantities such as the c entral 
black hole mass (ISeigar et al.ll2008t iBerrier et al.ll2013l) . 

It has long been known th at pitc h angle does depend on the 
distribution of mass (e.g.. ISeigar et al.il2006l. 120141) and on the 
size of the central bulge (for instance, the obse rved correlation 
bet ween pitch angle and sigma reported in lSeigar et al.l (120081) 
and IBerrier et al. l ( 2013 ). as well as the qualitative relation of 
pitch angle to bulge size featured in the Hubble classification). 
In addition, it has been reported that pitch angle var ies with 
the total mass of gas in galactic disks dRobertsIl 19751) . This 
Letter demonstrates the fundamental way in which we can un¬ 
derstand how the spiral structure depends on the distribution 
of mass in disk galaxies. Furthermore, it illustrates how the 
qualitative Hubble morphological types can exhibit varying 
pitch angles for galaxies that have similarly sized bulges and 

14 This is true for the total gas density: however, the atomic gas density 
generally has a peak value at some radius and decreases toward the center 
due to conversion of the atomic to molecular gas. 
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Figure 2. Three-dimensional plot of the plane defined by the fit of Equation 0 with the multivariate normally distributed sampling depicted by a translucent 
blue meshed surface, along with the plotted points of the 24 galaxy member included data set (depicted by red spheres with the Milky Way in green). Note 
that the points will appear slightly darker when they are projected behind and partially obscured by the plane. The axes [x,y,z] depict [tan |P|, log (M*“ ,sc /Mq), 
respectively. Left: the view has been oriented parallel to the plane. Middle: the view has been oriented at an orientation sufficient to view the face of the 
plane. Right: the view has been projected along an orthogonal vector above the plane. 


are thus categorized as the same type. These galaxies likely 
have different gas densities in their disks. Although the den¬ 
sity wave theory provided the inspiration for this study, other 
theories may also be able to explain this result. Certainly, any 
successful theory of spiral structure must be able to do so. 
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